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EXISTENCE AND REGULARITY OF POSITIVE 
SOLUTIONS TO ELLIPTIC EQUATIONS OF 
SCHRODINGER TYPE 



^ . B. J. JAYE, V. G. MAZ'YA, AND I. E. VERBITSKY 

H ' 

Abstract. We prove the existence of positive solutions with op- 
timal local regularity to the homogeneous equation of Schrodinger 
type, 

— div(^V7i) — (Tu = in il, 
under only a form boundedness assumption on a € I?' (51) and 
n , ' ellipticity assumption on ^ G L°°(r2)"^", for an arbitrary open 

set 17 C R". 

We demonstrate that there is a two way correspondence between 
. the form boundedness and the existence of positive solutions to this 

' equation, as well as weak solutions to the equation with quadratic 

nonlinearity in the gradient, 

~div{AVv) = (AVv) • Vu + CT infl. 

(N ; As a consequence, we obtain necessary and sufficient conditions for 

both the form-boundedness (with a sharp upper form bound) and 
the positivity of the quadratic form of the Schrodinger type oper- 
' ator Ti = — div(^V-) — a with arbitrary distributional potential 

. a G D'{Q), and give examples clarifying the relationship between 

' these two properties. 

o ■ 

1. Introduction 

^ ■ 1.1. The goal of this paper is to present an existence and regularity 

, theory for positive solutions to the equation of Schrodinger type: 

(1.1) - div(^Vu) - (TM = in fi, 

on an arbitrary open set Q C R", n > 1, under the standard ellip- 
ticity assumptions on ^ G L°°(f2)"'^", and the sole condition of form 
boundedness on the real- valued distributional potential a G D'{Q): 



(1.2) 



I (a, h^)\<C [ (AVh) ■ Vhdx, for all h G C^{n). 
Jn 



Date: January 13, 2013. 

Key words and phrases. Schrodinger equation, positive solutions, form bound- 
edness, elliptic regularity, weak reverse Holder inequality. 

The first and third authors are supported in part by NSF grant DMS-0901550. 

1 



2 



B. J. JAYE, V. G. MAZ'YA, AND I. E. VERBITSKY 



Simultaneously, a corresponding theory will be developed for (pos- 
sibly sign changing) weak solutions to the equation with quadratic 
growth in the gradient: 

(1.3) - dw{AVv) = (AVv) - Vv + ainQ. 

In displays fll.ip - (ll.3l) . A : Q ^ j^nxn ^ ^^^^ n x n (possibly non- 
symmetric) matrix-valued function on Q, so that there exist m, M > 
such that, for almost every x G 

(1.4) m|eP < A{x)^ ■ e, and \A{x)C\ < M\C\, for all ^ e R^ 

It has been a long standing problem to extend the existing theory 
to general classes of cr, including highly oscillating, singular or distri- 
butional potentials, where the separation of the positive and negative 
parts of a is impossible due to the interaction between them. In our 
framework of distributional cr, positive solutions to the Schrodinger 
equation are not locally bounded, and consequently standard PDE 
tools based on Harnack's inequality, and the classical iterative tech- 
niques of Moser jMosGOj and Trudinger |Tru73] , or their extension due 
to Brezis and Kato |B K79j . are no longer available. 

A primary result of the present paper is the following principle: For 
any form bounded potential a with the upper form bound strictly less 
than 1, one can find a positive solution of U.l\) which lies in the local 
Sobolev space This regularity for positive solutions of f ll.ip 

is in fact optimal in the generality of the potentials considered here, 
as demonstrated by examples discussed below. Furthermore, there is a 
two way correspondence between the existence of positive solutions u G 
of fll.ip satisfying an additional logarithmic Caccioppoli-type 
condition, and the form boundedness of the potential cr; see Theorem 

As a consequence, necessary and sufficient conditions will be estab- 
lished for both the form-boundedness and the positivity of the qua- 
dratic form of the Schrodinger type operator H = — div(^V-) — cr 
with arbitrary distributional potential a G D'{Q). The form bound- 
edness property (11.21) is known to be equivalent to the boundedness of 
the operator T-L : Lq ^(fi) — )■ L~^'^(f2) from the homogeneous Sobolev 
space Lq ^(fi) into its dual. It is therefore a natural class of poten- 
tials in which to study the Schrodinger equation. In a wide class of 
domains fl, it has been characterized by the second and third authors 
|MV02al[MV06] . 

Our results for the equations (II. ip and (I1.3P in turn provide an alter- 
native proof (with a sharp upper form bound) of the characterization 
of (ll.2p established in |MV02at IMV06] , where harmonic analysis and 
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potential theory methods were employed. In addition, we obtain a 
characterization of potentials a G T>'{VL) satisfying the corresponding 
semi-boundedness property, so that the operator "H is non- negative: 

(1.5) {a, h^) < [ (AVh) ■ Vhdx, for all h e C^in). 

Jn 

Both equations (11.11) and (11.31) . as well as the quadratic form prop- 
erties of "H (II. 2p and (II. 5p are of fundamental importance to partial 
differential equations, spectral theory, and mathematical physics. Con- 
sequently, these questions have attracted the attention of many au- 
thors, starting from the foundational work of Bocher, Hartman, Hille, 
and Wintner on the Sturm-Liouville theory (see e.g., |Hi48j . |Har82j . 
Chapter 11), followed by contributions of Agmon Ag83 , Aizenman 



and Simon [AS82j. Ancona |An86j . Brezis and Kato |BK79j . Chung 
and Zhao jCZ95] . Maz'ya [Mazllj . Murata |Mur86j et al. in the multi- 
dimensional case. A recent survey of this rich area has been given by 
Pinchover [PmOTj. We also refer to |Maz69l [CFKS891 INP92I IBNV94I 
IBM971 [RS981 IRSS94I ISW99l iFitOOl iMSOOl IShaOOl Mm02\ [DN021 [DD03] 
and references therein for equation (11.11) and form-boundedness prop- 
erties ([I2D, (O]), and |AHBV09l lE^ IFM98[ IFMOO] for equation 

Given the wealth of the previous literature, it is important to stress 
what is novel about our approach. In all the papers listed above, var- 
ious assumptions on the potential a ensure the validity of Harnack's 
inequality for positive solutions of the Schrodinger equation or some 
form of compactness properties of "H. Moreover, a is usually decom- 
posed into the sum of its positive and negative parts: cr = (T+ — (T_, 
which are treated separately, with more stringent assumptions on 0"+ 
than (T_ . In many of these results ct+ is assumed to belong to the Kato 
class of potentials, while ct_ to the local Kato class (see [AS82j . [CZ95j ). 
The corresponding positive solutions are continuous, and the existence 
of a positive solution is equivalent to the positivity of the quadratic 
form of H. In the mathematical physics literature the latter is known 



as the Allegretto-Piepenbrink theorem (see e.g. |CFKS89j . Sec. 2). All 



these tools are not available for general potentials a. 

The primary technical hurdles of our approach in comparison with 
the existing literature arise from the following essential characteristics 
of a satisfying (II. 2p : 

(1) a in general does not lie globally in a dual Sobolev space, i.e. 
a ^ L"^'''(fi) for any s > 0; 

(2) there are no local compactness conditions on a. 
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From the first item above, it is clear that one cannot follow standard 
methods to achieve global estimates which would yield the existence 
of solutions of (11. ip . Indeed, there are simple examples of cr so that a 
solution u of (11.11) does not lie in L^{Q). On the other hand, as a result 
of the second item, finding the correct quantity to work with in order 
to prove local estimates becomes a subtle issue. 

We will see that the two inequalities contained in (II. 2p are responsible 
for two distinct aspects of the existence of solutions to (11.11) and (II. 3p . 
Let us therefore consider the following upper and lower bounds of the 
quadratic form {a, h?): 

(1.6) ((J, h^) < \ I {AVh) ■ Vhdx, for all h e C^{n), 

Jn 

and 

(1.7) (a, h^) > -A / (AVh) ■ Vhdx, for all h e C^in). 

Jn 

In what follows, a positive function u on is defined to be a function 
u G such that u > quasi-everywhere in f2. Let us now state 

our first theorem: 

Theorem 1.1. LetQ C R", n > 1, be an open set. Let a G D'{Q), and 
A : Q j^nxn ^ matrix function satisfying the ellipticity conditions 
( [i.y^l ). Then the following statements hold: 

(i) Suppose cr obeys U.6\) with an upper form bound A < 1, and (1.1) 



with a lower form bound A > 0. Then a G L^^^'^iyL), and there exists a 
positive solution u G L{^^(VL) of the equation 

(1.8) - div(^VM) = au m D'{n), 

so that the following logarithmic Caccioppoli inequality holds: 

(1.9) [ ^-^^if'^dx <Co [ |V<^Pc/x, for all cp G Co"^(fi), 
Jn Jn 

with Co = Co{n, m, M, A) > 0. 

(ii) Suppose a G L^^'^(f2), and there exists a solution u G Lj^^(fi) 
satisfying ( fi.g|) for a constant Co > 0. Then there exists a solution 
V G L\'^^{yt) of the equation 

(1.10) -dw{AVv)=A{S7v)-Vv + a m D'{n), 
such that 

(1.11) [ \Vv\^ip'^dx < Co / \Vip\^dx, for all ^ G C^{n). 
Jn Jn 
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(iii) Suppose M.IC^) has a solution v G Lj^^(f2) satisfying U.ll\) for a 
positive constant Cq > 0. Then a satisfies the lower form hound ( [i. 
for a positive constant A = A(Co, m, M) > 0, and the upper form hound 
[rd\) with 

(a) X = 1 if A is symmetric; 

(h) \ = [ — if A is non- symmetric. 
\m J 

In the case where cr is a positive measure, the relationship between 
positive superharmonic supersolutions of f ll.Sp and the vahdity of f ll.6p 
has been discussed in a general framework using probabilistic meth- 
ods by Fitzsimmons in |FitOOj . This work was in turn building on 
the important paper of Ancona |An86] . In our framework, as we are 
considering oscillating potentials, one cannot rely on the theory of su- 
perharmonic functions, and we need to prove sharper estimates in order 

1 2 

to obtain the stronger L^^^ regularity. Without this it is not obvious 
how to even make sense of solutions to the equation ( II. ip . 

The sharpness of Theorem II. II is exhibited by well known examples. 
One such example is included in Sec. [7]for the convenience of the reader. 
Here it is also shown that in general there exist positive solutions u of 
(11. ip which lie in the space L\^^{Vl), but not L\^^{yt), and that statement 
(i) of Theorem 11.11 fails in general if A = 1. One should also note that in 
statement (iii), in the case of a non-symmetric matrix A, the constant 
A must in general depend on M and m, see Sec. [71 

1.2. In Theorem II. H it was seen that the lower form bound (I1.7P 
on the potential a is necessary in order to obtain solutions satisfying 
the regularity conditions (II. 9p and (II. lip . These conditions are of 
importance in our application of Theorem 11.11 in characterizing the 
inequality (II. 2p . and are also of classical interest in partial differential 
equations. 

However, if one is solely interested in the existence of solutions to 
(II. ip and (II. 3p . then the conditions on both the operator A and a may 
be relaxed to local conditions. We say that A satisfies a local ellipticity 
and houndedness assumption if, for each f/ CC fi, there exist positive 
constants mu and M(/, such that for any a; G f/, we have 

(1.12) mu\i\^ < A{x)^ ■ and \A{x)^\ < Mu\^\, for all ^ G R". 

The condition (I1.7P can be relaxed to a local condition stated in terms 
of dual Sobolev spaces: when n > 3, suppose 



(1.13) 



a = div(G) in n, where G G £i„""^((])". 



6 



B. J. JAYE, V. G. MAZ'YA, AND I. E. VERBITSKY 



Here Cf^^ ^(^2) is the local Morrey space defined in Sec. [2l see (12 ■2p 
and (12. lip below. This condition is significantly weaker than (II. 2p . In 
dimensions n = 1, 2 we only require a G L[~^'^(i7). 

The following theorem should be compared to statements (i) and (ii) 
of Theorem 11.11 above: 

Theorem 1.2. Let Q C R" be an open set, and suppose that A satisfies 
the local ellipticity and houndedness conditions U.1S\) . Let a G V{Q) 
satisfy U.6\) with a constant < A < 1, and in addition suppose 

(a) n = 1,2, and a G L^^'^(n); 

(b) n>3, and a satisfies the local condition 1^1. 13\) . 

Then there exists a positive solution u G L\j^^{yi) of M.l\) . and a solution 

veLlUn) ofm. 

A crude sufficient condition for (I1.13P is a G L|^(,^(fi) + Lj~).'"(f2), 
but much more general a are admissible for (ll.lSp . We emphasize that 
condition (I1.13P is considerably weaker than the usual local Kato class 
condition. It is not necessary for the existence of a positive solution 
u G of (II. ip . However, it is the sharp condition to obtain 

solutions u of (II. ip so that log(M) G BMOiod^)- 

To prove Theorems 11.11 and 11.21 we make crucial use of certain 
Caccioppoli-type inequalities. As was mentioned above, the classical it- 
erative techniques used in [MosGOl ITru73l [BK791 RS82l ICFG861 IMZ97j 
are not available. 

Instead, we interpolate between a Caccioppoli inequality and an es- 
timate on the mean oscillation of the logarithm to obtain uniform dou- 
bling properties on an approximating sequence. See Proposition 13. 8| 
which constitutes a key part of the argument. From this doubling prop- 
erty, one can deduce local uniform gradient estimates. This technique 
yields the optimal regularity for solutions of (II. ip in the generality of 
potentials satisfying (II. 2p or (11.130 . 

Along the way, we obtain a characterization of when a nonnegative 
weight function satisfying a weak reverse Holder inequality is doubling 
(see Sec. 2.2 for definitions). Our main hard analysis tool here is Propo- 
sition which may be of independent interest. 

1.3. Let us now turn to discussing applications of Theorem 11.11 As 
a first application, we deduce an alternative approach to the results 
of the second and third authors in |MV02aj . regarding the character- 
ization of the inequality (II. 2p . It avoids the heavy harmonic analysis 
and potential theory machinery that was used in |MV02aj . and IS con- 
siderably more elementary. This program is carried out in Sec. HI In 
particular, if = R" and with A the identity matrix, we will show 
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that the form boundedness condition (11.21) is equivalent to the follow- 
ing representation of a: 
(1.14) 

a = div(f), with [ h^\f\'^dx < Ci [ \Vh\'^dx for all h G C^{n). 

Moreover, ( 11.14p with Ci = \ implies (II. 2p with C = 1. Conversely, 
(II. 2p with C < 1 (or more precisely, (II. 6p with A < 1 and (I1.7P with 
A > 0) implies ffTTTip . 

In Section [5l we consider distributions a G V'{Q) satisfying the semi- 
boundedness property (II. 5p . In the case of the Laplacian, it means 

(1.15) {(r,h^)< [ \Vh\^dx, for all h e (Q) . 

Our main result in this regard is the following. (See Theorem 15 . 1 1 below 
for a similar criterion concerning the general operator div^(V-) in place 
of the Laplacian.) 

Theorem 1.3. A real-valued distribution a G V{Q) satisfies U.15\) if 
and only if there exists T G Lf^^^Q)"', so that 

(1.16) a < div(f) - |f|2 zn V'{n). 

The inequality in (I1.16P can not in general be strengthened to an 
equality. Such conditions have their roots in classical Sturm-Liouville 
theory in one dimension, see e.g. |Har82] (Chapter 11, Theorem 7.2). 

It had been conjectured that a condition characterizing (I1.15P was 
the following: 
(1.17) 

a < div($), where [ \h\'^\$\^dx <C I \Vh\'^dx, for all h G C^{n), 

for some $ G Lj^Q^(fi)" and C > 0. In other words, this means that 
'half the condition (11.140 found in |MV02aj should characterize semi- 
boundedness. However, it is proved below that, for any C > 0, condi- 
tion (I1.17P is not necessary for (11.150 to hold, although it is obviously 
sufficient when C = |. 

Proposition 1.4. Let Q = R", n > 1. Let a be the radial potential 
defined by 

n — 1 2 
a = cos r H sm r — sm r, 

r 

where r = \x\. Then a satisfies ^1.15\) . but cannot be represented in the 
form 
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This is the content of Proposition 17.11 in Sec. 7, where additional 
examples are exhibited to help clarify our results. 

Theorem 11.31 above concerned the case when A = 1 in (11. 6p . so 
that the Schrodinger operator fails to be coercive in the homogeneous 
Sobolev space. In the case when a is a positive measure, one can in- 
stead consider superharmonic supersolutions of the equation (II. ip in 
this critical case. This is a sharpening of |FitOO] mentioned above. 
Indeed: 

Proposition 1.5. Suppose that Q is an open set, and let a be a positive 
Borel measure defined on Q. Let A be a symmetric matrix function 
satisfying (T^. Then a satisfies 

(1.18) [ h^da < [ AiVh) ■ Vhdx, for all h e C^{n), 

Jq Jn 

if and only if there exists a positive superharmonic function u so that 

— div^(VM) > au in Vt. 

In Proposition II. 5[ the notion of superharmonicity is the one asso- 
ciated to the operator A, see e.g. |HKM06] . This is proved in Section 
ED below. 

In Section [6], we consider a recent resuh of Frazier, Nazarov and the 
third author |FNV10j on positive solutions with prescribed boundary 
values. Let us recall one of their main theorems. Suppose f2 is a 
bounded NTA domain, and that a is a nonnegative measure in f2. 
Then, under precise necessary and sufficient conditions on a up to 
the boundary, a positive minimal solution u (called the gauge in the 
probabilistic literature, see e.g. |CZ95] ) to the equation 

I —An = ou in f2, 
I n = 1 on oil, 

is constructed in |FNV10j (see Theorem 16.11 below). The solution is 
understood in the sense that 

u{x)= / G{x,y)u{y)da{y) + l, 
Jn 

where G{x, y) is Green's function of the Laplacian. In this paper, we 
will adapt the approach taken in the proof of Theorem II. II to show that 
in fact u G see Theorem 16.21 below. This regularity is again 

optimal under the assumptions of the theorem. 

In conclusion, we remark that our approach outlined above is non- 
linear in nature, and an extension to general quasilinear operators of 
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p-Laplacian type will be presented in a forthcoming paper |JMV12] . 
where the L^-analogue of f ll.2p will be characterized. 

1.4. Acknowledgement. We would like to thank Yehuda Pinchover 
for suggesting that Theorem 11.21 should require only local ellipticity 
and boundedness conditions on the operator A. 

2. Preliminaries 

2.1. Notation and function spaces. For an open set Q C R*^, n > 
1, we denote by C^{fl) the space of smooth functions with compact 
support in Q. The energy space Lq is then the completion of 
C^{Q) with respect to the Dirichlet norm || V/i||^2(-q). The majority 
of estimates in this paper are local; we say that h G if h(p G 

lI'^{Q) whenever (p G C^{Vl). 

For test function arguments it will be useful to introduce the space 
Lj'^(f2). We say that h G if /i G Lq'^{Q) has compact support. 

Define to be the dual of Then a distribution a G 

L;Jf{Q) if ifa G for any if G C^i^). 

We will write V CC U, for two open sets U,V C R", if there exists 
a compact set K C R" so that V C K C U. 

Throughout the paper, we use the usual notation for the integral 
average: 




For an open set U, we say u G BMO{U) if there is a positive constant 

Du so that 

(2.1) 

-f \u{y) — -f u{z) dz\^dy < Du, for any ball B{x, 2r) C U. 

JB{x,r) J B{x,r) 

In addition, u G BMOiodyt) if for each compactly supported open set 
U dd VL, there is a positive constant Djj > so that (12.11) holds. 

Let us next introduce the local Morrey space: we say / G £f^^(f2) if, 
for each compactly supported set U CC fl, there exists a constant Cu 
so that 

(2.2) [ \ffdx< Cus'', for all balls B{z, 2s) C U. 

Jb{z,s) 

We conclude with the definition of a multiplier (see |MSh09] ). Let X 
and Y be two normed function spaces, and let Z be a dense subset of X. 
We say that g is a multiplier from X to Y, written a.s g d M{X — )■ Y), 
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if (? ■ / G F for all / G and there is a positive constant C > so 
that the following inequality holds: 

\\9-f\\y<C\\f\\^, for all /GZ. 

In what follows X and Y will be Lq''^{^1) and respectively, and 

Z will be C^{n). 

2.2. On weak reverse Holder inequalities and BMO. In this sec- 
tion we characterize the weak reverse Holder weights that are doubling. 
This forms a key tool in our argument. First, let us introduce some 
notation. 

Definition 2.1. Let U C R*^ be an open set, and let w he a nonnega- 
tive measurable function. Then w is said to be doubling in U if there 
exists a constant Au > so that. 



(2.3) 



-+- w dx < Au4- wdx, for all balls i?(x,4r) C U. 

J B{x,2r) J B{x,r) 



Let w be a nonnegative measurable function. Then w is said to satisfy 
a weak reverse Holder inequality in U if there exists constants g > 1 
and Bjj > so that, 

(2.4) (i w'^dxY''< Bui wdx, for all balls B{x,2r) C U. 

Remark 2.2. The following simple consequence of the doubling prop- 
erty will prove useful. Let \J be an open set, and suppose w is doubling 
in f/. Then, whenever B{x, 4r) C U and z G B[x, r) with B{z, 4s) C U, 
we have 



w{y)dy < C{Aij, s,r)j- w{z)dz. 

B{x,r) J B{z,s) 

This principle will be used in a Harnack chain argument in Proposition 



Our argument hinges on the following result: 

Proposition 2.3. Let U be an open set, and suppose w satisfies the 
weak reverse Holder inequality ^2.4^ in U. Then w is doubling in U, 
i.e. holds, if and only if\og{w) G BMO{U) (see [2l\)). 

In particular, if w satisfies ( [^./^| ) and 

(2.5) f I \ogw{y) - f \ogw{z)dz\'^dy < Du, 

Jb(x,s) Jb{x,s) 



w dx, 
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for all balls B{x,2s) C U. Then there is a constant C{Bu,Du) > 0, 
so that for any ball B{x, 4r) C U 

(2.6) I wdx <C{Bu,Du)-f 

JB{x,2r) J B{x,r) 

where Bu is the constant from ^2.4\) - 

Only the the sufficiency direction is required in what follows; how- 
ever, since this characterization does not seem to appear explicitly in 
the literature we prove the full statement. To prove Proposition 12. 3^ 
we use the following lemma: 

Lemma 2.4. Let U C R" be an open set. Suppose that there exist 
s > 1 and w > 0, along with a constant Ci > so that the following 
inequality holds: 



w^dxj < Cij- wdx, whenever B{x,2r) C U. 

B{x,r) ^ J B{x,2r) 



Then, for any t > 0, there exists a constant Ct = C{t,Ci) > so that 

\i/s / r \ i/t 

dxj < Cti-h w^dxj , whenever B{x,2r) C U. 

B{x,r) ' ^JB{x,2r) ' 

This lemma had been used in proving estimates for quasilinear equa- 
tions by G. Mingione |Min07j . A proof can be found in Remark 6.12 
of |Giu03] . Let us now turn to proving the proposition. 

Proof of Proposition \2.3[ Let us first prove the necessity, suppose that 
w satisfies the weak reverse Holder inequality (12. 4p . and in addition 
that w is doubling in U. Then, for each ball B{x,Ar) C U, we have 

„ /• r 

w^dx\ < Bu+ wdx<AuBu4- wdx. 

B{x,r) ' JB{x,2r) J B{x,r) 

It follows that w satisfies a reverse Holder inequality in U, and is there- 
fore a Muckenhoupt Aoo-weight. It follows (see Chapter 5 of |St93] ) 
that log(M) e BMO{U). 

Let us now turn to the converse statement. Suppose w satisfies (12.51) 
and (12. 4p . From (12. 5p . it is a well known consequence of the John- 
Nirenberg inequality that there exists a constant < t < 1 so that 
is an 742-weight in U, i.e. there exists a positive constant A > 
(depending on Du in (12. 5p ) so that for all balls B{z, 2s) C U, 



(2.7) / w'dx<A(-f 

Jb(z,s) VS 



w * dx 

B(z,s) V_B(2,s) 
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Indeed, the John-Nirenberg inequahty (see |St93] ) yields t = t{Djj) > 
so that 



+ exp ( t 


\og{w){y') - f \og{w{y))dy 


Jb{z,s) ^ 


Jb(z,s) 



(2.8) 

The inequality (12. 8p contains two inequalities: 

/ exp(\og{w\y')) - / \og{w\y))dy)dy' < C{Du), and 

J B{z,s) ^ J B{z,s) ' 

I exp(log(w-*(2/')) + / \og{w\yy)dy\dy' < C{Du). 

J B{z,s) ^ Jb{z,s) ^ 

Multiplying these two inequalities together, one obtains ( 12. 7p . 

Combining (12. 7p with Jensen's inequality, we see that if -8(2;, 4s) C 
[/, then 



(2.9) / w'dx<Ard w-'dx) ^<A2"/ 

Jb(z.2s) Jb(z.s) ' JBi 



dx. 

< B{z,2s) J B{z,s) ' JB{z,s) 

Let B{z,8s) C U, then applying Lemma [2.41 yields 



-I wdx < Cu,t (/ dx) ^ < Cu,t (-[ dx 

Jb{z,2s) ' yB(z,As) ^ yB{z,s) 

< Cu.t-r wdx. 



'B{z,s) 

The second inequality in the chain follows from the doubling of w^, 
and the last inequality follows from Holder's inequality. By a standard 
covering argument, the factor of 8 in the enlargement of the ball can be 
replaced by 4, which yields (12.60 . This completes the proposition. □ 

2.3. Preliminaries for distributional potentials a. Let Q be an 

open set in R", with n > 1. Let A : Q ^ j^nxn ga^igfying the local 
conditions (11.121) . For a real- valued distribution a defined on fl, we 
define the multiplication operator by 

{ah, h) ■.= {(T,h^), for all h G C^{n). 

Suppose now a satisfies (II. 2p for a positive constant C > 0, and let us 
write: 

\\h\U= (^J (AVh) ■ Vhdxy^\ for h e C^{n) 

By polarization, we see that (11.21) is equivalent to the inequality: 

(2.10) \{crg,h)\ < C\\g\U\\h\U for all g,he C^{Q), 

with C > the same constant that appears in (11.20 . Furthermore, by 
the local boundedness assumption (I1.12p on the operator A, we can 
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extend f l2.10p by continuity, so that f l2.10p is valid for all g,h & Lq'^{U) 
whenever U GG Q. We denote this extension again by a. 

We now state some simple lemmas regarding the local character of 
the distributional potentials we will consider. Let us begin with an 
alternative way of stating the local condition f ll.lSp : 

Lemma 2.5. Suppose that n > 3 and a G V{Q). If a satisfies U.13\) 
then the following condition holds: 

For each open set U GG Q, there exists a constant Cu > 0, 

(2.11) such that whenever B{z, 2s) G U : 

\{h,a)\ < Cc7s("-2)/2||V/i|U2(B(,,,)), for all h G C^iBiz,s)). 

The two conditions f l2.1ip and f ll.l3p are in fact equivalent, but we 
omit the proof of the converse statement to Lemma 12. 5^ as we will not 
use it, and its proof is slightly lengthy. To prove the lemma, note that 
if 0" = div(G), with |G| satisfying fll.l3p . then it follows directly from 
differentiation that (12.1 ip holds. 

We next prove a key local property of distributions a satisfying (II. 2p . 
namely that (II. 2p is a stronger condition than the condition (12. lip . We 
define the capacity cap(-E', Q) of a compact set E by: 

(2.12) cap(E,fi) = inf{||V/i||i2(f^) : h G C^{n), h > 1 on E}. 

Lemma 2.6. Suppose that a satisfies and let U GG V GG Vt. 

Then a G L-^^^{U), and 

(2.13) ||a||L-i,2(c;) < MyC- cap(t/,l^)^/^ 

Here My is as in U.l^) . and C is the constant from U.^) . In partic- 
ular, for each ball B{x,r) so that B{x,2r) G U , we have 

(2.14) \ML-^a(Bi.,r)) < C^{C,Mu)r^^-^^'\ 

From display (I2.14p . it follows that (II. 2p is stronger than the local 
condition (EHH). 

Proof. Let h G C^{U), and let g G C^{V), so that ^ = 1 on [/. Then 
applying (I2.10p and (11.21) . we obtain 

\{a,h)\ = <C||(7|U||/i|U<MyC||V(7||2||V/i||2. 

Therefore a G L~^'^([/), and minimising over such g yields (I2.13P by 
definition of capacity. The second estimate is a special case of (12.130 
and follows from well known estimates for the capacity of a ball, see 
for example [Mazllj . □ 
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It will be convenient to use a mollification of the potential a. Let 
us fix a smooth radial approximate identity 9?, i.e. ip G C^(i?(0, 1)), 
so that (y9 > on -Bi(O), with Hv^Uli = 1- For e > 0, we denote 
ifi; = e~^(p{x/e). Then, denote the convolution of the distribution by 
cjg = v?£ * cr. The next two lemmas show that the mollification does not 
effect a in terms of form boundedness. We will write das = CTsdx. 

Lemma 2.7. Let Q be an open set, and let V CC Q. In addition 
let e < d{y,dVL) 12. Suppose that a G ViVt) so that M.(^) holds for a 
constant A > 0. Then the following inequality holds: 

(2.15) I h^dae < \ [ (AeVh) ■ Vhdx, for all h G C^iV). 

Jv Jv 

where Ae{x) = {ip^ * A){x). 

Proof. Let h G C^{V). We first note that by the interchange of molli- 
fication and the distribution (see Lemma 6.8 of [LLOlj ) 

{a,ipe*h'')= I ips{t){(T,h{--tY)dt. 

J B{0,e) 

By elementary geometry, h{ - — t) G C^{Q) for all t G B{0,e), and 
hence 

We, /i^) < A / ifeit) ( / {A{x)Vh{x - t)) ■ Vh{x - t)dx)dt 

(2 IQ) JB{0,e) ^Jn ^ 

Aeix)Vh{x) ■ Vh{x)dx, 

which proves the lemma. □ 

Our second mollification lemma says that if a satisfies the local con- 
dition fl2.1ip . then so does the mollification of a. Let us introduce the 
notation 11^ = {x e W : dist(x, U) < e}. 

Lemma 2.8. Suppose n > 3 and a G is a real valued distribution 

satisfying \2.11\) . Let V CC ^, then if e < d{V,dVL/2), the mollified 
potential satisfies /12.11\) for all open sets U CC V , with constant 

Proof. Let U CC be a compactly supported open set. Then, for 
B{x,2r) CU andhe C^{B{x,r)), h{-~x) G C^{n) for all x G S^(0). 
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Hence, for all such h, it follows that 
hda, 



e 

n 



< I ips{x)\{a,h{- -x))\dx 



<C^,r("-2)/2||v/,||^,. 
This completes the proof. □ 

3. The proofs of Theorems 11.11 and 11.21 

In this section we prove our primary existence theorems, as well as 
prove the connections between the solutions of the equations (II. ip and 
(11.31) with the validity of (11.21) . Let n > 1. Throughout this section 
we will assume without loss of generality Q C R" is a connected open 
set; note that in the case of an arbitrary open set, our arguments apply 
to each connected component. This assumption is used in a Harnack 
chain argument. 

The most substantial argument will be the assertion of statement (i) 
in Theorem II. H and its local variant in Theorem 11.21 We restate these 
two results as propositions for convenience. In light of Lemma 12. 5| the 
existence result for equation (II. ip in Theorem 11.21 follows from: 

Proposition 3.1. Let A : Q j^nxn ^ possibly non-symetric ma- 
trix function satisfying ( fi.i^) . Suppose that a is a real-valued distri- 
bution satisfying the local dual Sobolev condition \2.11\) . and the upper 
boundedness condition M.6\) for a constant < A < 1 . Then there is a 
positive solution u G Lj^^(i7) of M.l\) . 

The second proposition concerns the case when A satisfies global 
ellipticity and boundedness, and a in addition satisfies (II. 7p : 

Proposition 3.2. Suppose A : ^ j^nxn. ^ possibly non-symetric 
matrix function satisfying jl-j^ - In addition, suppose that a is a real- 
valued distribution satisfying ( [i. 7| ) for a positive constant A > 0, and 
M.0\) for a constant < A < 1. Then there is a positive solution 
u G L\'^^{yi) of M.l\) . Furthermore, u satisfies 

(3.1) I ^^<^2 dx<C [ I V<^P dx for all cp G C^(n). 

Jn Jn 

It was seen in Lemma 12.61 that if a satisfies (I1.6P and (II. 7p . then 
the condition (12. lip holds. Hence the existence part of Proposition 
follows from Proposition 13.11 
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3.1. An approximating sequence. To prove Propositions 13. ll and 
13. 2^ we use local properties of a to find solutions to a mollified variant 
of equation (II. ip in a sequence of subdomains of VL. We will then prove 
a uniform gradient estimate on this sequence. 

Suppose that a satisfies (11. 6p with < A < 1. Let for j > 1 be 
an exhaustion of VL by smooth domains, that is, Vtj GG ^j+i, and Q = 
[jj flj. In addition, let us fix a ball B so that its concentric enlargement 
AB CC Qi. Let Eq = 1, and Ej = mm{Ej^i/2,d{Qj,dQj+i)/2,2~^) for 
j > 1. With this notation, define: 

(Tj = (fsj * c", and Aj = (fg. * A. 

Referring to (I1.12p . we let rrij = ttiq. and Mj = Mq. be the local 
elhpticity and boundedness constants relative to Qj. For k > j, note 
that Ek < d{Qj, dflj+i), and hence for any x G Qj we have 

(3.2) m,+i|^p < Akix)^ ■ and \Akix)^\ < M,+i|e|, for all ^ G R'^. 

Define Uj to be the solution of 

\ —div{ A jVuj) = ajUj in Qj, 

u^dx = 1, Uj > q.e. 

Furthermore, Uj satisfies the Harnack inequality in Qj. 

Proof of existence and uniqueness of liS. 3\) . We will see that the exis- 
tence of (13. 3p is a simple consequence of the Lax-Milgram lemma. De- 
fine a bilinear form C on Lq'^^Qj) x Ll''^{Qj) by 

C{w, h) = {AjVw) ■ Vhdx — {(JjUJ, h), with w,h E Lq ^(fij). 

By the assumptions on a, the following properties hold: 

\C{w,h)\ < (M,-+i + C(n)||(T,|Uoo(j,^,)mf/")||Vw;||2l|V/i||2, and 

C{w, w) > mj+i(l — A) II Vw||2 . 

The first inequality follows from (13.20 along with the Sobolev inequality 
(this is true for all n > 1, using standard Sobolev inequalities). The 
second inequality is a combination of Lemma [2.71 and (13. 2p . Hence the 
hypotheses of the Lax-Milgram lemma are satisfied. 

Applying the Lax-Milgram lemma, we see that there exists a unique 
Wj G Lq'^IQj) satisfying 

C{wj, h) = {aj, h). 
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Let Vj = Wj + 1. Let us next show that vj > q.e. To see this, let 
(f = Vj, and note that (p~ = min(t>j,0) G LQ'^(f2j). By testing (13.31) 
with the vahd test function (f~, it follows: 

(3.4) / AjiVif') ■ Vf' = {ajif~,if') < A / Aj{Vip') ■ Vf'. 

Since < A < 1, it follows from (11.41) that (p~ = q.e as required. Let 
us now define 

Uj={j-^Vjdx^ Vj. 

Then Uj solves (13. 3p . The validity of Harnack's inequality for Uj follows 
from classical elliptic regularity theory, see e.g. |Tru73] . since o"-, is 
smooth. □ 

3.2. Caccioppoli and Morrey estimates for the approximating 
sequence. We next turn to proving two estimates on the gradient 
of the approximating sequence. The first estimate is a Caccioppoli 
inequality: 

Lemma 3.3. Suppose that a satisfies U.6]) with < A < 1, and let 
{uj} be the sequence constructed in ^3. 3\) . Let ip G C^{flj), then for 
any k > j , 

|2 „/,2 ^ r^/n^ ™ \\ I „,2 



(3.5) / \Vukr^lj'dx<C{Mj+i,mj+i,X) / ui\V^\'dx. 

Proof. Let us fix k and j as in the statement of the lemma, and let 
V = Uk. With ip G C^{Qj), > 0, test the weak formulation of (13.31) 
with vip'^ G Lq ^(i7j). Using (II. 6p . it follows that 



{{AkVv) ■Vv)^/j^ dx < {o-kV,ip^v) + 2Mj+i / v^p \Vv\\Vi{j\ dx 
= {ak{ipv),iljv) + 2Mj+i / vi/j\Vv\\Vip\dx 

<X Ak(y{vij))(y{vij))dx + 2Mj+i vij\Vv\\Vij\dx. 

Here we have used Lemma [2.71 and (akVjip'^v) = {o'k{vip),vip) . Using 
the Cauchy inequality, it follows that for any e > there exists a 
constant Cs, depending on e, A, Mj+i and rrij+i, such that 



(1-A) / {AkVv)-Vv^'dx<e \Vv\'^p' dx + v'\Vilj\'dx. 
Choosing e < {1 — A)mj_|_i and rearranging, we recover (13. 5p . □ 
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The second estimate we use relates a bound on the gradient of the 
logarithm independent of j, with uniform properties on the negative 
part of the quadratic form associated with the distribution a: 

Lemma 3.4. With {uj)j the sequence constructed in jiS. 3\) . the follow- 
ing estimate holds for all k > j and ip G C'^{VLj), 



(3.6) 



< -C(Mj+i,mj+i) I ^^dak 
+ C{Mj+i,mj+i) / iVtpl'^dx. 



Proof. Let h = ip'^/uk, with ip G C^{Qj), > 0. Since satisfies the 
Harnack inequality in Qj, there exists a constant c > so that Uk > c 
on the support oi ip. It follows that h is a valid test function for the 
weak formulation of fl3.3p . This yields 

(3.7) (AkVuk) -v^^yx = -(^ak,^^y 

On the other hand, 

m,+i [ ip'^dx <— [ (AkVuk) ■ v(—^dx 

+ 2M,+i / \Y^\V^\^dx, 

and therefore Cauchy's inequality yields a constant C = C(M,_|_i, rrij^i) > 
such that 



\Vuk\' 



'j^dx<-C [ {AkVukyv(—)dx + i 



(3.8) / ^ — ^ip'dx<-C I {AkVuk}V{ — )dx + C I \V^\'dx 
Combining (13. 7p and (13. Sp . we deduce (13. 6p . □ 



From Lemma 13. 4[ we deduce two estimates, depending on the addi- 
tion properties of a. First, we deduce that the so-called logarithmic 
Cacciopolli inequality holds if a in addition satisfies (II. 7p . 

Lemma 3.5. Suppose that the real-valued distribution a satisfies ji. 7| ) 
with constant A > 0, and let {uj)j be as in i^3. 3\} . Let ip G C^{VLj), 
then for any k > j 

(3.9) [ ^^^T/^^rfx < C(M,+i,m,+i,A) /" \V^/j\^dx 

Jo ui Jo 
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Proof. From Lemma 13.41 it clearly suffices to show that, for all ip G 
C^{flj) we have 

But this follows in precisely the same manner as Lemma 12.71 □ 

Recalling the definition of capacity in f l2.12p , we arrive at the follow- 
ing corollary of Lemma 13.51 

Corollary 3.6. Suppose the hypotheses of Lemma \3.5\ are satisfied. 
Then there exists a positive constant C = C{A,mj^i, Mj+i), so that 
whenever F CC Qj, we have 

(3.10) [ ^-^^dx < Ccap(F, Qj+i), for all k > j. 

Jf 

In the case when a only satisfies fl2.1ip . a local Morrey space estimate 
holds, which is a weakened version of (13. 9p : 

Lemma 3.7. Suppose n > 3, and that a satisfies \2.11\) . Consider the 
sequence {uj} as in hS. Then for each ball B{x, r) so that B{x, 2r) C 
flj, it follows that for all k > j 

|2 



(3.11) / < ^(^.^ mj+i)r"-2 



Proof. Fix such a ball B{x,r) as in the statement of the lemma, and 
let ip e C^{B{x,2r)) such that = 1 on B{x,r), < ip < 1 and 
{Vipl < C/R. The lemma follows from estimating (13. 6p with this choice 
of ip. For A; > j, it suffices to prove that 



(3.12) / ip'dak <C{nj)r 

Picking U = Qj in the definition of (12. lip , Lemma 12.81 yields 



ip'^dai, 



The last inequality follows since < ip < 1. Here {^j)ej is the ej- 
neighbourhood of Qj. Note that {^j)ej C ^j+i, and since Ek < Sj, 
it follows {^j)ek C {^j)ej and hence by definition in (12. lip we have 

The display (13.120 now follows from the estimate on the gradient of 

i). □ 
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In the case = 1 or 2; note that if a G ^[^^'^(fi), then for all k > j 
we have 
(3.13) 

f \^ k\ ^ C(fij, Mj_|_i, m-j+i) whenever B{x,2r) C Qj. 

J B{x,r) '^k 

In fact this estimate holds for all dimensions, but it is not strong enough 
to provide us with a uniform bound in higher dimensions. The estimate 
( ]3.13p follows from display (13. 6 p in Lemma 13.41 Indeed, for k > j, 
just pick the test function ip in (13. 6p so that ip = 1 on Qj, and ip ^ 
Co~(^+i). This yields 

r \Vuk\ ^ C(^.,Mj+i,mj+i, ||cT||i-i.2(o^,_^^)). 

Here we are using the fact that the mollification does not effect the 
local dual Sobolev norm within Qj, which can be established precisely 
as in Lemma [2^81 The estimate (I3.13P clearly follows from the previous 
display. 

3.3. A local gradient estimate. The key estimate is the following: 

Proposition 3.8. Suppose a is a real-valued distribution defined on 
Q, satisfying ( fi.61) with < A < 1, and in addition suppose that ^2.11\) 
holds. Let {uj} he the sequence in Then, whenever u CC VLj, 

the following estimate holds: 

(3.14) / \Vuk\^dx < C{u, A, A, mj+i, M,+i, B, Qj), for all k > j. 

J UJ 

Note that the estimate (13.141) is independent of k for k > j. This 
is the key to allow us to deduce the existence of positive solutions to 

m- 

Proof. Fix j, and k, as in the statement of the proposition, and let 
V = Uk- It suffices to prove that whenever B{x, 8r) CC Qj, there exists 
a positive constant C > 0, depending on n, A, A, rrij+i, Mj+i, B, 
B{x,r) and Qj such that 



(3.15) / 

JB(x'i 



\Vv\'^dx < C. 

{x,r) 

The reader should keep in mind that all constants will be independent 
of k. Fix such a ball B{x,8r) GG Qj. To prove (13.151) . we will employ 
Proposition 12.31 in U = Qj to show that v"^ is doubling in Qj, with 
constants independent of k. To verify the hypothesis of Proposition 
we first show that f ^ satisfies a weak reverse Holder inequality, i.e. 



EXISTENCE AND REGULARITY OF POSITIVE SOLUTIONS 21 

the inequality (12 .4^ holds in Qj. To this end, let us fix B{z, 2s) CC Qj. 
Suppose first n > 3. Let ip G C^{flj), then applying the Sobolev 
inequality yields 

(3.16) ( / v^\ij\^dx)~<C \Vv\^ij^dx + C v^\Vij\^dx. 

Applying Lemma 13.31 in the first term on the right hand side of fl3.16p 
results in 



(3.17) ([ v^\ij\^dx) " <C [ v^\Vilj\'^dx. 

We now specialise fl3.17p to the case ip ^ C^{B{z, 2s)), with ip = 1 in 
B{z, s), and < C/s. As a result, we obtain 



n-2 

(3.18) (/ {v^)^dx) " <C-f v^dx. 

^Jb{z,s) ' Jb{z,2s) 

The constant in C > in f l3.18p depends on n, M^+i, mj+i, and A. 
Hence, if n > 3, (12. 4p holds in U = Qj, with w = v"^ and q = n/{n — 2). 

If n = 2, we slightly modify the above argument. The following 
Sobolev inequality is standard (see e.g. |MZ97j . Corollary 1.57): for 
each g < oo, and for all / G C^{B{z, 2s)), 

(3.19) (/ \f{y)Vdy)"''<C{q){f \V f {y)\'' dy)"\ 

^JB{z,2s) ^ ^Jb{z,2s) ^ 

Using f l3.19p as in (13.161) and following the argument through display 
f l3.18p . it follows in the case n = 2 that (12. 4 p holds in U = Qj, with 
w = v"^ for any choice g < oo. Note that in the case n = 1 even 
stronger Sobolev inequalities are at our disposal, and so the estimate 
(12. 4p continues to hold; we leave this to the reader. 

To apply Proposition 12. 3[ it remains to show log(f) G BMO{Qj). 
For this let B{z,2s) C Qj. The Poincare inequality yields a constant 
C = C{n) such that 

\\ogv--f \ogv\'^ dx < Cs^-"" [ ^-^^dx. 



'B{z,s) JB{z,s) JB{z,s) % 

First suppose n > 3. Then from Lemma [3.71 we have 

|2 



(3.20) / < C{Mj+i, rrij+i, fijOs"-^ 

Jb{z,s) % 



and hence, 



(3.21) i \\ogv-4- \ogv\'^dx <C{Mj+i,mj+i,nj) 

Jb(z,s) Jb{z,s) 
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In the case n = 1,2; we apply the weaker estimate f l3.13p in com- 
bination with Poincare's inequahty, to conclude that fl3.2ip remains 
true in these cases. From (13.211) . we conclude (see (12.11) ) that \ogv G 
BMO{Qj), with SMO-norm depending only on n,mj^i, Mj^i,Qj. In 
particular, v"^ satisfies both (12.41) and (12. 5p in Qj. From Proposition 
it follows that v'^ is doubling in flj, with constants depending on 



n, m 



1, Mj+i, A and A, see (12. 6p . 



Since ilj is a smooth connected set, one can find a Harnack chain 
from B{x,2r) to the fixed ball B CC fli. In other words, there are 
positive constants Cq, Ci and > 0, depending on the smooth param- 
eterization of flj, along with points Xq, . . . x^- and balls B{xi, 4rj) C flj 
so that 

(1) B{xo,ro) = B{x,2r), and B^xnjTn) = B; 

(2) ri > Comin(ro,rAr), and |5(xi,ri)n5(xi+i,ri+i)| > Ci min(ro, tat) 
for alH = . . . A^ - 1. 

Since v"^ is doubling in Qj, a Harnack chain argument (see Remark 12.20 
applied to the chain construction above yields 

+ v'^dx < C{B{x,r),mj^i, Mj^i,ilj, B, X, A)-!- v^dx. 

JB{x,2r) Jb 

It therefore follows from the normalization on f ^ that 

(3.22) -I v'^dx < C{B{x, r), m^+i, M^+i, ilj, B, A, A). 

JB{x,2r) 

To complete the proof, combine the Caccioppoli inequality (Lemma 
13. 3p with the estimate (I3.22p . This results in the inequality 

C f 

Wvl^dx < — / v^dx < C, 



2 

B(x,r) ^ J B{x,2r) 

for a constant C > 0, depending on n, rrij^i, M,_|_i, B, A, A, VLj and 
B{x,r). Hence (13.151) is proved for a constant independent of k. □ 

3.4. Proof of Propositions r3.1l and l3.2i We begin by proving Propo- 
sition |3lTl before moving on to prove Proposition [3]2l and with it state- 
ment (i) of Theorem 11.11 

Proof of Proposition \3.1[ Let Qj be an exhaustion of Q by smooth do- 
mains. We will use Proposition 13.81 repeatedly in each Qj to deduce 
the existence of a solution of (11.10 . Fix a ball B G Qi, with AB C Qi, 
and note that the construction of the approximate sequence from (13.30 
is valid under the present assumptions on a, and the gradient estimate 
flXIil) holds. 
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First, by ( 13.141) with j = 1, along with weak compactness and Rel- 
hch's theorem, we pass to a subsequence u^^'' of uj so that u^p — )■ u^^^ 
weakly in Llf^{ili), and u^^^-' — )■ u'^^'^ both in L^(fii) and almost every- 
where in Let ej^i be the corresponding sequence from (13.31) . Since 
cr G it follows that whenever h G C^(fii) 

(3.23) {cruf\h) = {a,ufh) {a,u^^^h) = {au^^\h). 

Note also, by combining the uniform bound (I3.14p . with convergence 
of the mollification, we have 

\{a,^,uf\h) - {aufM < \\oe,,,-o\\L-.ai^^^^\\V{ufh)\\L2 

< C\\cxe^ -, - o-| 1^-1,2(^1) as j oo. 

We conclude: 

(3.24) {(Xj^iuf \h) ^ {au'^^\h). 

Similarly, by linearity and local boundedness of the operator A, we 
deduce that 

AVu^P -Vhdx ^ I AVu^^^ -Vhdx 



From the uniform bound (I3.14p . along with standard convergence prop- 
erties of the mollification in L^: 



— 7- as j — > cxD. 



/ {Aj^i- A)Vuf -Vhdx 

It therefore follows that the limit function u^-^^ G satisfies 

(3.25) - div(^VMW) = an(i) in ©'(^i). 

Given {u\^^}k, let us apply estimate (I3.14p in Qj+i to obtain a subse- 
quence u\^~^^'' of M^f^ and u'^^^^^ G Ll^l{Qj+i) with: 

^k^^^ u^^^^^ weakly in Ll^l{Qj+i), and u'^^^^^ u'^^^^^ a.e. in Qj+i. 

Note that by Lemma 12. 6[ it follows that a G L~^'^{Qj). As in the 
argument leading to (13.250 . we see that u^^~^^^ satisfies 

(3.26) - div(^VM(^+i)) = aw^^+i) in P'(^],+l). 

By construction, u^^^ = u^^'^^^ in Qj. Hence one can define u G 
by: u = u^^^ in ilj. By (13. 26^ it follows that 

-div(^VM) = au in V'{n). 
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Next, let us demonstrate that u is not the zero function. To see this 
note that 

(3.27) j {uffdx = 1, for all j, i. 

Since uf^ — > u in Li^^{Qi), we may pass to the limit in f l3.27p . A 
standard application of Mazur's lemma shows that the limit solution 
u > 0. On the other hand, for any k > 0, it follows from Lemma 13.41 
and weak compactness that there exists v G so that \og{u^^''^) — )■ 

V almost everywhere. But then v = log(M) a.e. and therefore (see 
for example Theorem 1.32 of |HKM06] ) log(^^^) converges weakly to 

log('u) in Ll^^{Qk)- Hence u > quasi-e very where, and m is a positive 
weak solution of (11. ip . □ 

We now move onto Proposition 13.21 Recall that here the matrix A 
satisfies the global ellipticity and boundedness conditions, so in all the 
previous estimates of this section, we replace rnj = m and Mj = M. 

Proof of Proposition \3.S[ Let us keep the notation from the proof of 
Proposition 13.11 The existence of a positive solution u G LI^^{^1) of 
(II. ip follows from Proposition 13.11 It was proved above in addition 
that log(M) is well defined in and in each log(M) is the 

weak limit of a sequence logw^. From Lemma [3.5^ it follows that, for 



k\2 



(3.28) / -^^tjj'dx <C{M,m,A) / dx. 

Since Vlogwj^^ converges to Vlogw weakly in Ll^^{Vtk), we deduce from 
weak lower semi-continuity of the norm that for all ip G C^{VLk) 

(3.29) /" ^^^^^rfa; < C(M,m,A) ! \Vil)\^ dx. 
Jn Jn 

For any k > 0, the estimate (13.290 holds with a uniform constant for 
smooth functions supported in Q^, and hence (13. ip holds. □ 

3.5. A logarithmic change of variable: solutions of ( 11.31) from 
solutions of ( II. Ij) . This section is concerned with deducing solutions 
of (II. 3p from solutions of (II. ip by a logarithmic substitution. This 
substitution is classical, for instance it appears in the study of ODEs 
in |Hi48] . and there are examples that show it can be delicate, see e.g. 
|FMOO] . In |AHB V09] . there is a rather comprehensive account of the 
connection between these two types of equations when cr is a finite 



EXISTENCE AND REGULARITY OF POSITIVE SOLUTIONS 25 

measure. We will prove the following lemma, from which statement 
(ii) in Theorem 11.11 and the remainder of Theorem 11.21 follows. 

Lemma 3.9. Let Q be an open set, and let A : Q ^ j^nxn gg^Ugjy 
the local conditions 1^1.12) . Suppose that a G L^^'^(fi), and that there 
exists a positive solution u of M.l\) . Then v = log(M) G is a 

solution of U.3\) . 

Proof. Fix U GG Q. The first step is to prove that 

(3.30) / ^-^dx <C(U,a,n,p). 

Ju « 

Let e > 0, and let V be such that U GG V GG n. For h G C^{V), 
test the weak formulation of (11. ip with ip = h{u + e)~^ G Ll''^{V). This 
yields 

(3.31) / ^ L.\/hdx= I — '——hdx+{a ,h). 



u + e Jn {u + u + e' 

Let us now estimate the third term on the right. By assumption a G 
L~^''^{V), and hence there exists f G (L^(\/))" so that a = div(r) in 
V. Therefore 

(3.32) {a^^,h) = [ ^^^(^—]hdx+ f ^^\/h-fdx. 

u + e u + e \u + eJ J^u + e 

Since 6/{u + e) < 1, it follows from Cauchy's inequality that for any 

5 > 0, we have 

(3.33) 

\{a^^,h)\<6 f -J^^hdx + Cs f \f\^hdx+ f \Vh\\f\dx. 
u + £ Jq{u + er Jq Jn 

Now, let G C^{V), > 0, = 1 on f/, and put h = in flCTl) . 
Rearranging, using the ellipticity and boundedness assumptions (II. 4p . 
we obtain 
(3.34) 

|2 



iTT'V I 1 —Lp^dx<2Mv / \V'^\u)dx 



2 



+ ^ / T dx + Cs / \T\^ip^dx + 2 / |V<^| |r|(^cia; 



n 



Here the bound fl3.33p has also been used. Appealing to Cauchy's 
inequality again in (I3.34p . we obtain 

|2 
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Letting e — )■ 0, f l3.30p follows from Fatou's lemma. 

Now, let us again look at fl3.3ip . this time with an arbitrary h G 
C^{U). It follows from (ICTj) that as £ 



nu + e Jn u 

(AVu) ■ Vm , , f (AVu) ■ Vn , , 
— — — nax ^ / h ax. 

To handle the last term in f l3.3ip . note that from f l3.30p and the domi- 
nated convergence theorem 



'( — )= ( ^ in L\^) as £ ^ 0, 

It 



on the other hand, it is clear that — )■ 1 in L (U)^ as £ — ?■ 0. Thus, 

u 

-i+e 

conclude that 



it follows that -)■ 1 in L^''^{U) as e ^ 0. But since a e L ^'^{U) we 



(a , h) = (a, h) — )■ (ex, h), as £ — )■ 0. 

u + e u + e 

It follows that V = log(M) is a distributional solution of (11. 3p . □ 

Proof of Theorem \l.l\ statement (ii). This is nothing more than a re- 
statement of Lemma 13.91 above, along with the trivial observation that 
if u satisfies (11. 9p . then v = log(M) satisfies (II. lip . □ 

We may now also complete the proof of Theorem II. 2[ 

Proof of Theorem By Lemma 12.51 and Proposition 13. 1| it follows 
that under the hypothesis of Theorem 11.21 there exists a position so- 
lution u G of (II. ip . By Lemma [3^ setting v = log('u), we see 
that V G is a solution of (II. 3p . □ 

We did not use (ll.7p or global ellipticity and boundedness assump- 
tions in the previous lemma, doing so allows us to conclude that the 
solution satisfies an additional multiplier condition. 

Lemma 3.10. Under the assumptions of Lemma \3.y[ suppose that A 
satisfies (T^, and a satisfies ( [i. 7p for a positive constant A > 0. Then 



there exists a solution v G of l{1.3\) satisfying 

Proof. We will keep the notation from the proof in Lemma 13.91 It is 
left to prove that v satisfies (II. lip . To this end, let us again test (II. ip 
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with '^'^ l[u + e\ for a smooth test function Lp G > 0. There 

exist constants Ci, C2 > depending on m and M, so that 

(3.35) 

(m + e)^ \ V + ^ V + 

Indeed, as in display fl3.3ip . one obtains by testing equation fll.ip the 
following identity: 



Jn u + e ^ Jn {u + ey ^ ^ u + e'^' 
Hence, from the ellipticity and boundedness assumptions (11. 4p . we have 

f IV?/P f IV?/I V 

(3.37) m/ .^^L\ydx<M .\V^\dx - {a^,^'). 

J^{u + eY Jnu + e u + e 

From an elementary application of Cauchy's inequality in (I3.37p . dis- 
play fl3.35p follows. One can pick, for instance, Ci = (M/m)^ and 
C2 = 2/m. 

Next, applying (11.70 . it follows that the second term on the right 
hand side in (13.350 is bounded by a constant multiple of 

(3.38) A / \Y^(^)\^dx + A [ ^\Vip\'dx. 

\u + eJ ^ Jnu + e^ 

The first term in (I3.38P converges to zero as e — ?■ 0, by virtue of (I3.30p 
and the dominated convergence theorem. Again by dominated conver- 
gence, the second term in (I3.38P converges to \ V(p\'^dx, as e — )■ 0. 
Substituting these estimates into display (I3.35p . we deduce that (II. lip 
holds. □ 

3.6. Existence of solutions to ( II. 3j) implies the validity of ( II. 2|) . 

Proof of Theorem statement (iii) . This will follow immediately from 
Lemmas 13.111 through 13.131 below. □ 

Lemma 3.11. Let Q be an open set, and suppose that a is a distri- 
bution defined on Q. Let A be an n x n real-valued symmetric ma- 
trix defined on Q satisfying ([i.^[ ). Suppose there exists a supersolution 
V G L]'^^{9) of mE), then (EJ) holds with A = 1. 

Proof. Suppose that there exists a solution v of (11.30 . and let ip G 
C^(fi), then testing ([L3D with yields 

(ayp, yp) = (a, < 2 / \^\\{A{\/ v)) ■ \/ ^\dx - [ (AVv) ■ {Vv)<p'^dx. 

Jn Jn 
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Under the present assumptions, ^ is a symmetric positive definite ma- 
trix. It follows that for C,,i] & R", 

(3.39) \iAO-v\<iiAO-0'^'iiAv)-vy^'- 

Thus, we have that 

(3.40) 

2 [ \ip\\{AiVv))-Vifi\dx<2 [ \ip\{{AVv)-VvY/^{iAVip)-VipY/^dx 
Jn Jn 

< [ ip^iAVv) ■Vv)+ [ {(AVip) ■ Vip)dx. 
Jn Jn 

Therefore dH]) holds with A = 1. □ 

On the other hand, if symmetry is not assumed, one may still con- 
clude the validity of (11. 6p and (I1.7P from (II. 3p . as we will show now. 

Lemma 3.12. Let Q be a connected open set, and suppose that A : Q ^ 
R"^", satisfying Suppose that v G L]f^{Q) is a supersolution of 

mM), then [rd\) holds with 

V m / 

Proof. Let us first show that (II. 6p holds with the given choice of A. To 
this end, let ip G C^{VL) and test the weak formulation (II. ip with the 
valid test function ip^. Together with the assumptions (II. 4p . this yields 



(crv9, v?) = (cT, v?^) < M2 / \V v\\lp\\V Lp\dx - m I \Vv\^ip^dx 

Jn Jn 



(3.41) 



< 
m 



where Young's inequality was used in the last line. □ 



We conclude this section by showing that if the solution v of (II. 3p 
in addition satisfies (II. lip , then a satisfies (II. 7p for a positive constant 
A > 0. 

Lemma 3.13. Under the assumptions of either Lemma \3 . 1 1\ or Lemma 
\3.12l if one in addition assumes the solution v G L]^1{Q) satisfies 
U.ll\) . then a satisfies ( [i. 7| ) for a positive constant A > 0. 
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Proof. Let ip G C^{Q) and test the weak formulation (11. ip with the 
vahd test function Then 



for a suitable choice of constant A > 0. Appealing to (ll.4p . the result 



3.7. On the equation ( II. 3p . In this section, we make a few comments 
regarding our results for the equation (II. 3p in comparison to existing 
literature, as the introduction of this paper focussed rather more on 
the Schrodinger type equation f ll.ll) . First, we will restate a theorem 
which has been proved for reference: 

Theorem 3.14. Let Q be an open set. Let a G D'[Q) be a real-valued 
distribution. Suppose A : Q ^ R" is a symmetric real-valued matrix 
function satisfying The the following two statements hold: 

(i) Suppose that a satisfies ( [i. 7p for a positive constant A > 0, and 
M.(j^) for a constant < A < 1, then there exists a positive solution 
V G L]^1{Q) of ( |j.3|) so that ( li.i ij) holds. In addition, the constructed 
solution has the exponential integrability property: e'" G L]^1{Q). 

(a) Conversely, if there exists a solution v G L]^{Q) of U.3\) so that 
U.ll\) holds, then a satisfies ([i. 7p for a positive constant A > 0, and 
/ fO) with A = 1. 

In |FM98j . solutions of (11.30 are proved in the global energy space 
LQ'^(f2) are proved when is a bounded domain, under the assumption 
that a G They explicitly note that this condition on a is used 

to guarantee that (ll.2p is valid. Theorem 13.141 therefore compliments 
their theorem with a more local result in nature, and therefore one 
which requires less restriction on a. As was noticed in [FMQSj IFMOOj . 
there exist classes of solutions of (11.31) that are exponentially integrable. 
One can trace this principle back to the employment of certain non- 
linear test functions in proving the existence of solutions to (ll.3p (see 
e.g. |Ev90l IFMOOj ). A refinement of this argument is what is also 
employed in the current paper, since we deduce Theorem 13.141 from 
our considerations of the Schrodinger type equation via a logarithmic 
substitution. 

The local exponential integrability in statement (i) is sharp, as can 
be seen from the example discussed in Section [71 The paper [FMOO] 



(3.42) 




follows. 



□ 
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concerns quasilinear equations of p-Laplacian type, and we will consider 
such equations in our forthcoming paper [JMV12j . 

3.8. On the critical case A = 1. In this subsection we briefly discuss 
the limiting case when A = 1. We shall prove the following proposition: 

Proposition 3.15. Suppose that Q is an open set, and suppose A is 
a symmetric matrix function satisfying (T^- Then U.18\) holds if and 
only if there exists a positive superharmonic function u such that 

(3.43) - div(^Vn) > au in f2. 

Proof. Let us assume that VL is connected. We shall assume n > 2, 
and leave the one dimensional case to the reader. The necessity is well 
known, and holds even in very general potential theoretic frameworks, 
see e.g. [Fit 00]. For the converse, let \j G (0,1) be a sequence such 
that \j — 7- 1. Applying a very special case of Theorem 11.11 above, we 
find a sequence of positive functions {uj}j G L\j^^{yi) of 

(3.44) — div(^VMj) = XjO'Uj, with + Uj dx = 1. 

Jb 

Here B CC is a fixed ball. We may assume that Uj is superharmonic. 
Next, let us fix a smooth connected subdomain U of Q, so that B CC U. 
For a fixed q < n/{n — l), we will prove that for any ball B{x, r) CC U, 
the following estimate holds 

(3.45) / \Vuj\'^dx <C{m,M,U,q,B{x,r)). 

JB{x,r) 

To this end, note by the property (11.91) that 

(3.46) [ /i^ dx < C{m, M) [ \Vh\'^dx for all h G C^i^). 

Display f l3.46p is also a standard property of superharmonic functions, 
as a consequence of Moser's work |Mos60] . From f l3.46p and the Poincare 
inequality, one readily deduces as in Proposition l3.8l that Uj G BMO{U). 
Therefore, as in Proposition 12.31 from the John-Nirenberg inequality 
we find a constant c = c{m,M,n) > 0, < c < 1, such that u'^ is 
doubhng in U (see (12. 3p ). with doubling constants depending on m, M, 
and n. For any B{x,r) CC U, a Harnack chain argument yields a 
constant C = C{{B{x,r), U, B, m, M) such that 

(3.47) f u^dx <C I u'dx < C{B{x, r), U, B, m, M), 

JB(x,r) Jb 

where in the last equation, we have used normalization of Uj in f l3.44p . 
Let us now note that following inequality superharmonic functions. 
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essentially due to Moser |Mos60] : for 1 < q < n/{n — 1), there exists a 
constant C = C{n,q) such that 

(3.48) [ iVujl" dx < Cr""-" inf uj. 

JB{x,r) -S{a;,r) 

Combining f l3.48p and fl3.47p . the display f l3.45p follows. Note that from 
f l3.48p and Rellich's theorem, we deduce that there exists u such that 
Uj — > u a.e in fi, and j^^udx = 1. Furthermore, u can be chosen to be 
superharmonic, by standard convergence properties, see |KM92] . 

Our aim is to show that u = lim inf j.^. q.e. To this end, let 
V = liminij^oo Uj, and denote by v* the lower semi-continuous reg- 
ularization of v. By the fundamental convergence theorem for super- 
harmonic functions, v* is superharmonic, and v* = lim inf j_5.oo uj quasi- 
everywhere (see Theorem 7.4 of |HKM06] ). Since v* = u a.e. and they 
are both superharmonic, we have that u = v* everywhere. The claim 
follows. 

Let us now conclude the argument. Since a satisfies fll.l8p . it does 
not charge sets of capacity zero. By Fatou's lemma, for any G C^{VL) 
with > it follows that 

liminf / Ujipda> / u^pda. 
Jn Jn 

Combining this last display with the weak convergence of the Riesz 
measure for superharmonic functions, see e.g. |HKM06] . we conclude 
that ([S3SD holds. □ 

4. Form boundedness 

In this section we apply Theorem 11.11 to deduce a new proof of The- 
orem 14.11 below, which was the primary theorem in |MV02a] (see also 
[MV06j ). Let us consider the case Q = IV^, n > 3, since the case of a 
general domain, under certain mild restrictions on Q, can be reduced 
to the entire space, as was explained in |MV02aj . 

Theorem 4.1. Let a G i^'(R"), n > 3. Then the following statements 
hold. 

(i) The quadratic form inequality 

(4.1) h')\ < C ||V/i||i.(K„), for all h G Co^(R"), 
is valid if and only if a can be represented in the form 

(4.2) cT = divf. 
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where T G ^^^^(R")" obeys 

(4.3) / /i2|f|2rfx < Ci ||V/i||i2(R„), for allheC^iR''). 

(ii) If a satisfies fH?T]) . then (g^D holds with f = V(A"V), where 
A-V G BMO(R"), and 

(4.4) 

/ /i2|V(A-V)prfx<Ci||V/i||i2(R„), for allheC^iK'). 

(iii) // gSD holds with Ci = I then dH]) holds with C = 1. Con- 
versely, if a satisfies (11. 6p wzt/i t^e upper form hound A < 1, and (11.71) 
with the lower form hound A > 0, then (14. 4p /ioWs wzi/i a constant Ci 
which does not depend on a. 



Proof of Theorem ^.i, The sufficiency part of statement (i) with Ci = 
C^/4 in inequahty (14. 3 p follows using integration by parts and Cauchy's 
inequality: if a = div F, then 



r ■ Vh h dx 



(4.5) 

To deduce the remainder of the theorem, we apply part (ii) of The- 
orem [LT] with a := o"/(A + e), where A is the upper form bound of a, 
and £ > 0, so that the corresponding upper form bound A of o" satisfies 
A < 1. This yields the existence of a weak solution \E' G L]'^^(R") of 
the multi-dimensional Riccati equation: 

(4.6) -A^ = |V*P + 5- in L)'(R"). 

Using /i^, where h G C^(R"), as a test function in this equation, and 
integrating by parts, we estimate 



h'\sm'dx = 2 -Whhdx- {a,h') 

R" Jr" 

< 2 I I W^l U2(R„) I I V/l| U2(Kn) + A I I V/l| |i2(R„), 

where A = A/(A + e) is the lower form bound of a. From this it follows: 

(4.7) f \V^!\^dx < (1 + V^)2 ||V/i||i2(R„), for all h G C^{YC). 
Jr" 

In other words, |V^| G M(Lo'^(R") ^ ^^(R")), where M(Lo'^(R") ^ 
L^(R")) is the space of pointwise multipliers from the homogeneous 
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Sobolev space Lq ^(R") into L^(R") defined in Sec. 12.11 Hence, a can 
be represented in the form 

a = -div - |V^|^ |V^| e M(Lo'^(R") ^ ^^(R")). 

Moreover, one can deduce from (14. ip that a is well defined in the 
sense of the weak-* BMO convergence (see details in |MV06j ): 

A-\^Ijn <?) A-^a e BMO(R"), 

where iPn{x) = ip{\x\/N), and ip G C^(R) is a standard cut-off func- 
tion. It follows from (g^D that A^^ (IV^H G BMO(R"), and hence 

(4.8) ^ = -A-^ (|V^n - A^V G BMO(R"). 
Thus, a can be represented in the form 

(4.9) a = divf, f = VA^^ a, in D'(R"), 
where 

(4.10) A"^ a = - A-^ (|V^n G BMO(R''). 

To complete the proof of statements (ii) and (iii), it remains to verify 

VA-^a G M(4''(R") ^ L2(R")). 

Let g = (— A)^3 |V\E'p > 0. In other words, g is the Riesz potential 
of order 1 of |V^'p, so that jVA'^ (|V^n(x)| < c{n) g{x) for almost 
every x G R". 

Since |V ^| G M(lJ'^(R") ^ ^^(R")), it follows (see Theorem 1.7 
in |Ver99] ) that g G M(L^'2(R") ^ ^^(R")). Hence, VA"^ (IV^H G 
M(Lo'^(R") ^ L2(R")). Thus, 

f = VA-^a G M(LJ'^(R") ^ L2(R")). 

This is equivalent to (14. 4p . and the proof is complete. □ 

5. Semi-boundedness 

Let fl C R" be an open set with n > 1, and suppose ^ is a matrix 
function satisfying ( II. 4p . In this section, we will consider real- valued 
distributions a G which are semi-bounded; that is, the quadratic 

form of the operator H = — div(^V-) — cr is non-negative: 

(5.1) {a, h^) < I {AVh) ■ Vhdx, for all h G C^in). 

Jn 
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It had been conjectured that a necessary and sufficient condition for 
(15.11) to hold is the following condition: there exist T G Lf^^{Q)"- and a 
constant C > so that cr < div(r), and 

(5.2) [ \h\^\f\^dx <C [ {AVh)-Vhdx for all h e C^in). 
Jn Jn 

A simple estimate using integration by parts and Cauchy's inequality 
in the form (I3.39P shows that condition (15.21) with C = | is sufficient 
for (15. ip to hold. However, it is not necessary, with any C > 0, for 
(15. ip . We defer the proof of this fact to Proposition 17.11 below. On 
the other hand, the following theorem provides a characterization of 
semi-bounded distributions. 

Theorem 5.1. Let Q be an open set, and let a G be a real valued 

distribution. In addition, let A be a symmetric matrix function defined 
on Q satisfying ( [i.^| ). Then i\5.1\) holds if and only if there exists a 
vector field T G L'^^^{Q) so that 

(5.3) a < div(^f) - (^f) ■ f zn V'{n). 

There is an extension of Theorem 15.11 to non-symmetric matrices 
A. Indeed the necessity of the condition (15. 3 p extends to the non- 
symmetric case, see Proposition 15.31 below. On the other hand, a rep- 
etition of the proof of Lemma 13.121 shows: if (15. 3p holds, then (15. ip 
holds with a constant (M/m)^ introduced in the right hand side. Here 
m and M are the ellipticity constants from (II. 4p . 

Remark 5.2. The proof of Theorem 15.11 shows that (15. ip holds if and 
only if there exist solutions to the differential inequality: 

(5.4) - div(^VM) - (AVu) ■Vu>a in Q. 

The inequality in (15. 4p cannot be strengthened to an equality for gen- 
eral distributions a G V'{Q). Indeed, if there exists a solution v G 
of the equation 

-div^Vn = (AVu) -Vu + a in 

then it follows that a G Lj~^'^(i7) -|- Ll^^{^l). For instance, when n > 3, 
one can pick a = —S^q for xq E fl, then obviously (15. ip holds but a does 
not lie in the aforementioned class. In fact, in the special case when a 
is a measure, it is known that a G L[^^'^(r2) + Ll^^{Q.) if and only if a 
does not charge sets of capacity zero (see Theorem 2.1 of |BG096j ). 

Let us now move onto proving the Theorem: 
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Proof of Theorem \5.1[ The sufficiency of f l5.3p for (15. ip is a repetition 
of the proof of Lemma 13. Ill The necessity of (15. 3p is somewhat more 
involved, and follows from Proposition 15.31 below. □ 

Proposition 5.3. Let Vt he an open set, with A a (possibly non- 
symmetric) matrix function defined on satisfying Let o G 
©'(fi) satisfying / (5. Then there exists T G Lf^^{Q)"' so that 

(5.5) a < div(^f) - (^f) ■ f zn V'{n). 

Proof of Proposition 15. 5t Without loss of generality, we may assume 
that Q is connected. Otherwise, we simply repeat the argument which 
follows in each component. 

Step 1 {Approximation). Let f2j be an exhaustion of f2 by bounded 
smooth connected domains. Let \j G (0, 1) be any sequence so that 
Aj — 1, and 

Ej < l/2min(d(nj,anj+i),2-^). 

Consider cxj = Xj^Pe^ * cr, and denote by Aj = ipe^ * A. Then from 
Lemma 12. 7[ it follows that (Tj satisfies 

(5.6) [ \h\^daj<\j [ Aj{Vh)-Vhdx, ioi ail h e C^{Vtj). 

In addition to ( 15. 6p . note the following estimate: for any h G C^{Qj) 

(5.7) / \h\'daj>~C\\a,\\L^(^n,)\n,n^h\\l, 

this is a consequence of standard Sobolev inequalities (for all n > 1). 
By (11. 4p . we conclude that 

(5.8) / \h\^dcxj > -C\\aj\\L^(^n,)\^j\^^"' I AjVh-Vhdx. 

From (15. 6 p and (15. 8p . we see that the hypothesis of Theorem 11.11 are 
satisfied in Qj with potential aj. It therefore follows from Theorem ll.il 
that there exists Vj G L^^^{Qj) so that 

(5.9) - dw{AjVvj) = AjVvj ■ Vvj + aj, in V'^Qj). 

By addition of a suitable constant, we may assume that, for a fixed 
ball B CC Qi, 



(5.10) 



Vjdx 



1, for all j. 



Step 2 {A uniform bound). Fix 1 < j < k. Our aim is to show that 
Vk G Ll^l{Qj), with constants independent of k. Let h G C^{Qj), by 
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testing the weak formulation of Vk in (15.91) with h"^, we deduce from 
n that 

m \Vvk\'^h^dx < M 2h\Vvk\\Vh\- h'^dak- 
J^k J^k 

Applying Cauchy's inequality in the first term on the right hand side: 

ml \Vvk\^h^dx <— f \Vvk\'^h^dx + 2— [ \Vh\'^dx- f h^dak, 
Juk 2 J^k ^ J^k J^k 

and hence, as \k G (0, 1), we have 

(5.11) / \Vvk\^h^dx<C [ \Vh\^dx + C\{ipek*h\a)\. 



Next from standard distribution theory (see e.g. |Str03j . Chapter 8), 
it follows that 

\{^,,*h\a)\<C, 

for a constant C depending on a, the support of h and | l^'^'^/il |loo for 
some collection of multi-indices ai, . . . ,a]sf- (One can see this from 
either the structure theorem, or by the definition of continuity). In 
conclusion, for any h G C^{flj) 

(5.12) / \Vvk\^h'^dx <C{a,h). 

Juk 

This proves the claim that Vk G Lj^l{^lj), with constants independent 
of fc. 

Step 3 {Conclusion). This will be quite similar to Section l374l In- 
deed, consider first Vti. Then from (15.121) and weak compactness, we 
find a subsequence Vj^i of Wj, and G L\^^{yii) so that Vj^i — >■ weakly 
in L[Qj,(f2i). From (15.101) and an application of Rellich's theorem, the 
limit function vi is not identically infinite. 

Let G G (L^(i7i))" n with compact support in Vti. Note that 
(5.13) 

/ Aj^iWvj^- G dx = i {{Aj^i- A)Vvj^i) ■ G dx + i AVvj^i- Gdx. 
Jn ' Jn ' ' Jn ' 

The first term on the right of (I5.13P converges to zero as A; — > oo. 
Indeed, one can estimate: 



/ iiAj,i-A)Vvj,i)-Gdx 
Jn 



<||G'||oo(/ \Vvj,i\^dx) ( \Aj,i-A\^dx^ 

supp(G) ' ^ ^isupp(G) 
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and the right hand side converges to zero by (15.121) and standard prop- 
erties of approximate identities. For the second term on the right hand 
side of (I5.13p . note that by weak convergence 



/ AVvj^i -Gdx ^ I AVv^ ■ G dx, 
Jn ' Jn 



as J — )■ oo. 



It follows that, for all G G (L^(fii))" fl L°°, with compact support in 

(5.14) / Aj^Vvj^i ■ Gdx I AVv^ ■ Gdx, as j oo. 
Jq ' Jn 

It is not difficult to see that one can extend (15.141) for all G G (L^(fii))", 
with compact support in Qi. In other words, that Aj^iVvj^i — t- AVv^ 
weakly in L^^^{Qi)"-. 

Let h G C^(f2i) so that h > 0. We next claim that 

(5.15) liminf / {Aj iVvj i) ■ Vvj ihdx> / (AVv^) ■ Vv^h dx. 

in ' ' ' Jn 

To see this, denote by A^ the symmetric part of A, i.e. 2A'^ = A + A^. 
Then as before it follows {Aj^iYVvj^i — )■ A'^Wv^ weakly in L'^^^iVti). In 
addition, by standard properties of mollification, {Aj^iY — )• A'^ in the 
weak-Tir topology of L°°(f2i). We will repeatedly use the observation 
that the non-symmetric part does not contribute toward the quadratic 
form. First note that as a result of the weak convergence, we have 



(5.16) 



f {AVv^)-Vv^hdx= [ (A'Vv^) -Vv^hdx 

= liminf / {Aj lYVvj i ■ Vv^hdx. 



Next, by symmetry of the matrix (see (I3.39P ). we estimate for each j 

1/2 



/ (^j, i)'Vfj, 1 ■ Vv^hdx <( {Aj^iY^Vj^i ■ Vvj^ihdx 

[ {Aj,iY^v^ -Vv^hdx 



1/2 



By taking the limit infimum of both sides, using the weak-T^r convergence 
of the convolution, one obtains f l5.15p . Here we are using the following 
elementary fact for two bounded sequences (a^) and (bj): 

If liminf aj = a > and lim bj = b > 0, then liminf ajbj < ab. 
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On the other hand, by standard properties of the convolution, and since 
A, ^ 1: 

Aj(cr, ipi.. * h) {a, h), as j — )■ oo. 
Keeping ( 15.9p and f l5.14p in mind, we conclude that 

(5.17) - diY{AVv^) > (AVv^) ■ Vv^ + a in ©'(Qi). 

For k > 1, and given the sequence {vj^k-i}, a repetition of the above 
argument yields a subsequence Vj^k of so that Vj^k converges to 

(5.18) - div(^Vi;^) > {AVv'') ■Vv'' + a in V'iQk)- 

Furthermore, as in Section [331 "we may assert that = v^~^ in Qk-i- 
One can therefore define a function v G so that 

-div(^Vt;) - (^Vf) -Vvya in V'{Q). 
To complete the proof, it suffices to let F = —Vv. □ 

6. The local regularity of solutions to the Schrodinger 
equation with prescribed boundary values 

The goal of this section is to apply the regularity techniques devel- 
oped in this paper to the recent work of Frazier, Nazarov and Verbitsky 
[FNVlOj . The point here is to prove regularity of a given solution with 
prescribed boundary values, when we already know there exists a ma- 
jorant of the given solution. 

Let Q C R" be a bounded domain so that the boundary Har- 
nack inequality is valid (for instance a Lipschitz, or more generally 
a NTA domain). Let us now fix xq G Q, and let G{x,Xq) be the 
Green's function for the Laplace operator relative to Q. Then we define 
m{x) = min(l, G(x, a;o)). If is a C^'^ domain then m{x) is pointwise 
comparable to dist(a;,9fi). 

Let 0" be a locally finite Borel measure in Q. Then the following 
theorem is proved: 

Theorem 6.1. [FNVlOj Suppose that a satisfies the following embed- 
ding inequality 

(6.1) [ h^da <\[ \Vh\^dx, for all h e C^{n), 
Jn Jn 

with < A < 1. In addition, suppose that there is a constant c > so 
that 

(6.2) / m(x)exp( — -— ■ / m{y) da {y)\ da [x) < oo. 
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Then there is a solution ui of the equation 

-Aui = aui in Q, 



(6.3) 



Ml = 1 on dQ. 



Conversely, if there is a solution of Ii6.3\) . then li6.S\) holds for a positive 
constant c = c{Q), and i\6.1\) holds with A = 1. 

The solution constructed in Theorem 16. II is interpreted in the poten- 
tial theoretic sense, i.e., Ui G L^{Q,mda), and 

ui{x) = / G{x,y)ui{y)da{y) + 1. 
Jn 

If is a bounded C^'^ domain then ui G L^{Q,dx) fl L^{Q,mda), and 
is a solution to (16.31) in the very weak sense (see [FNVIO] ). 
Our primary result in this section is the following: 

Theorem 6.2. Consider the solution ui of constructed in Theo- 
rem \6.1[ then 

(6.4) u, G L\;^S^). 

Proof. Let VLj be a exhaustion of Vt by smooth domains. Let ipj G 
C^{VL) be such that ipj = 1 on VLj and < < 1. Note that (as in 
Lemma [2^ ipja G L~^''^{VL). Since (16.11) holds with constant < A < 1, 
we follow a similar argument to Section 13. ![ using the Lax-Milgram 
lemma to obtain a unique Vj G L^'^(n) satisfying: 

— Atij = fjCrVj G Q, 

Note here that Vj-^i > vj. 

Let B{x,2r) CC Q. By repeating the proof of Lemma [373| we deduce 
that there is a constant C = C{n) so that 

(6.6) [ \^Vj\^dx<^! v]dx. 

JB{x,r) ^ JB{x,2r) 

From (16. 6p . one asserts, as in displays (I3.16P through (I3.18p . that 

(6.7) (/ vf/^-'Uxf~"^"<ci v]dx. 

Since (16. 7p holds for all balls -B(x, 2r) CC fi, we see that the hypothesis 
of Lemma [2.41 are valid. Applying the lemma with t < n/(?T, — 2), one 
finds a constant C(t) > so that 



(6.5) 



(6.8) (/ v]dx)'' <C{t){l v]dx^'' 

yB{x,r) ' yB{x,2r) 
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Combining (16. 8p with (16. 6p . we conclude 

(6.9) [ \Vv/dx<C{t)r''-'^(-f v'jdxY^'. 

Jb{x,t) V_B{a;,2r) ^ 

We now wish to show that Vj < Ui. For a locally finite measure a, 
denote by Q'^{x, y) the minimal Green's function of —A— a (see |FV10j ). 
i.e. the minimal positive solution u{-,y) of the equation 

-Au{-,y) - au{-,y) = 6y, y e Vt. 

Then Ml = 1 + e;-(x, y)da{y), and Vj = l + G^^'^ix, y)^j{y)da{y) 
(recall ipja G W^^'^ (Q), so this representation coincides with the 
unique solution). By construction Q'^ is monotone in a (it can be 
represented by a Neumann series), and since ipja < a, it therefore fol- 
lows that Vj < Ui, for each j. Here Ui is as defined in (16.31) . Thus, from 

(6.10) [ iVvA'^dx <Cr''^H—^ — ^ [ u{dxY^\ 

JB{x,r) ^|5(a;,2r)| 7B(x,2r) ^ 

Letting t < n/(n — 2), and recalling the weak Harnack inequality (valid 
since a > 0), we deduce the estimate: 

(6.11) f |Vt;jprfx < Cr"-2f inf ui)^ 

JB(x,r) \B{x,r) / 

Using (16. lip , we readily deduce that there exists v G L\j^^{VL), so that 

1 2 

Vj increases to f , and Vj — v weakly in L^^^. Furthermore, as in the 
proof of Proposition 13.21 

(6.12) -Av = av inV'iyt). 

Since Vj < mi, it follows v <ui. On the other hand, we see that v = 1 
on dVl in the potential theoretic sense. Indeed, with G'(yu) denoting the 
Green's function of the Laplacian relative to Vt applied to yU, we see by 
f l6:T2|) that 

(6.13) v = G{av) + h, 

where h is the greatest harmonic minorant of v. Clearly, h > 1. On 
the other hand, since f < Mi, 

h-l<v-l<ui-l = G{aui). 

Thus /i — 1 is a nonnegative harmonic minorant of G{aui). However, 
by the Riesz decomposition theorem, the greatest harmonic minorant 
of G{aui) is zero. Thus h = 1, and w = 1 on dfl in the potential 
theoretic sense. 
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By minimality of mi, it thus follows that v = ui, and ui G 



□ 



7. Examples 



Our first result in this section completes our discussion of the condi- 
tion (15. 2p . and that it does not provide a characterization of distribu- 
tions satisfying (15. ip . 

Proposition 7.1. Let Q = R", and A be the n x n identity matrix. 
Let a be the radial potential defined by: 

a{x) = COST -\ sinr — sin^r, r = a; G R"'. 



Then a satisfies ( 15. but cannot be represented in the form l^5.S\) . 

Proof. We first consider the case n = 1 and Q = R+ = (0, -|-oo). Note 
that a criterion of form boundedness takes the form |MV02aj . |MV02bj : 
a = T' where 



r 



(7.1) 




a > 0. 



Let a = cosx — sin^ x. Then a is semibounded by (15.31) . i.e. 




but a cannot be represented in the form 
(7.2) a = r' - /i, /i > 0, 



where F satisfies (17. ip with any C > 0. 



In fact, even a weaker condition 



(7.3) 




as 



a — )■ +00 



cannot be satisfied if a is of the form (17.21) 
Indeed, suppose 

(7.4) cos X — sin^ x = T' — fj,, 



/x>0. 



Then 



r(a;) = sin X (^1+2 
where ip{x) is nondecreasing on R_|_. 
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V5-1 



.904 < 1 so that cos a — sin^a > for 



Let ao = arccos ■ ^ 
— ao < a < ao, and consequently F is nondecreasing in the interval 
[27cn — ao, 27m + ao]. Hence, for a = it follows that 

r(a + 2?™) - r(27rn) 

1 



sin a I 1 H — cos a h (^{a + 27m) — ip{27m) 



> sin a ( 1 H — cos a\ = C 



where C is independent of n. Here C > since sin a; (l + | cos x) — 
is increasing on (— ao, ao) and equals zero at the origin. 
On the other hand, for a = a + 27m we have 



Via) < 



ao 



a+a 2 

r(a;) dx < — 
ao 

r'2-ITn 



ct+1 



\Tix)\ dx. 



Similarly, r(27rr2) > '^q j^^2Tvn^^^^ Hence 

2 

r(a+27rn)-r(27rn) < 



a+27m+l 



a+2iTn+l 2 

|r(a;)| dx^ — 

(^0 J -a+2-Kn 



\T{x) \ dx. 



(^0 J a+2-Kn 

By (17.31) the right-hand side of the preceding inequality tends to zero 
as n — )■ +00. This contradicts the estimate 

r(a + 27m) - r(27rn) > C > 

obtained above. 

We now are in a position to consider the multi- dimensional case 
= R", n > 3. Let 



(7.5) 



a = cos r + 



n — 1 



■ sm r — sm r, 



\x\ 



xeW' 



Then by (15.31) with F = ^ sinr it follows that (15. ip holds. 
Note that satisfies the inequality 



(7.6) 



2 smr 
h dx 



R" 



Indeed, using polar coordinates and integration by parts, we obtain 



Sin 7^ 

h{rO- r^'-^drd^ 



R4 



(7.7) 



< 



cosr 



2h{r^) Vh{r^) ■ ^ r"-^ drd^ 



R4 



+ {n-2) 



/ / cosr /i^(r^) r" ^ drdC, 
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Applying Cauchy's inequality and Hardy's inequality, we estimate the 
last line by 

2(f ^dxf\\Vh\\,.^^.) + {n-2)[ ^dx 

< C ||V/l||^2(R,n), 

and hence (17. 6p holds. 

It remains to show that a\ = cos r — sin^ r does not satisfy the 
inequality Ui < divF, for any F which obeys fl5.2l) . This is equivalent 
(see |MV02a] ). to finding F with divF is form bounded. We note 
that as was proved in |MV02a] . we may always pick F = V$ where 
divF = A$, and 



(7. 



8) / /i2|V$pdx<C||V/i||i2(R„), forall/iGCo"^(R" 



for some C independent of h. Suppose now that 
(7.9) (Ti < A$, 

where ^ satisfies (I7.8p . Since o"i is radially symmetric, it follows by 
using the average of $ over the unit sphere that (17.91) holds with a 
radially symmetric ^o{x) = so that 

7X — 1 

(^i{r) < (fi" (r) -\ ip'{r), r > 0. 

r 

Moreover, 



and hence 



From this and (17.81) with a radially symmetric test function h, we de- 
duce 

[ /i^(r) |^'(r)|2r"'idr < — J— /" [ \V<^ir^)\^ h^{r) r'^-^drd^ 

<C [ \h'{r)\^r''-Ur. 

Consequently, 

pa 

|(/?'(r)|2r"-ic/r < Ca"-^ a > 0. 
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We now let ilj{r) = (f'{r). It follows from the preceding estimate that 

fa+l ^ 

[r) \ dr < C a~2^ a > 0. 

J a 

This implies 

^^c/r<Ca-2, a>0. 
r 

Next, denote by g{r) the function: 

"+00 



g{r) = ip{r) — (n — 1) / dr. 



r 

It is easy to see that g satisfies the same inequality as ip: 



/ \g{r)\dr < C a'^ , a > 0. 

J a 



Furthermore, 



Tl — 1 Ti — 1 

g'{r) = tp'lr) H iplr) = ip"{r) H (p{r) = A$(x). 

Thus, cri(r) = cosr — sin^r < g'{r), r > 0, where g satisfies the 
condition 

ra+l 

\g{r)\ dr = o{l), a — )■ +00. 



J a 



Hence, by the one-dimensional example on R_|_ considered above with 
g = r, we arrive at a contradiction. □ 

In the following example we consider non-symmetric operators, with 
the aim to show that the non-symmetric part can effect the constant 
appearing in the form bound. In Theorem II. H it was shown that when 
one has a solution to the Schrodinger equation (II. ip . or Riccati equation 
(11.31) . whose operator A is non-symmetric, then a satisfies (II. 6p with A 
depending on the ellipticity constants (II. 4p . This example shows that 
such a conclusion is not artificial: 

Example 7.2. Let n = 3, and suppose A = I -\- B, where B has 
zero entries except for 61^2 = Ca{xi) for a constant C, and a{xi) a 
Lipschitz continuous function. Suppose in addition 62,1 = ~^i,2- If 
u{x) = 1 -|- Ixp, then u solves 

1 / N —6 + 2x2Ca'(xi) 

— div(^VM) = (7U{x), with a = 

1 -I- |xp 

On the other hand 

^(o-e = ier- 
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It follows that in the case of non-symmetric matrices A, the constant in 
statement (iii) of Theorem 11.11 depends on the constant C, and hence 
the operator A. 

The next example (which is well known) demonstrates the sharpness 
of our primary theorem for Schrodinger type equations. In particular 
we confirm the assertions made in the introduction. 



Example 7.3. Consider positive solutions u of 
(7.10) 



-Au = r^u in R", 
infRn u = 0, 



with c < (n — 2)^/4. It is well known that (17.101) has positive solutions 
(up to constant multiple) of the form u±{x) = la;]"*, where 



(7.11) «± = ^±^V(^-2)2-4c. 

If c < (n — 2)^/4, then by Hardy's inequality it follows that (II. 6p holds 
with < A < 1. For < c < (n — 2)^/4 we see that by choosing a+, 
there exists a solution m+ G of (I7.10p . Taking c arbitrarily close 

to {n — 2)^/4, we see that the existence of a solution n+ G L^^^{Q) of 
(17.1 op is the optimal local regularity. The same example shows that 
solutions need not be locally bounded, and therefore positive solutions 
of (I7.10p do not satisfy the Harnack inequality. 

Choosing a_ in (17. lip , it follows that there exist positive solutions 
M_ G Lioc(R'') of dno]), which are not in Lj^^(R"). 

Finally, let c = {n — 2)^/4 in (I7.10p . The resulting unique positive 
solution does not lie in LJj^^(R"). This latter point shows the assump- 
tion < A < 1 in (II. 6p is necessary in order to prove statement (i) of 
Theorem 11.11 We remark that the uniqueness of the positive superhar- 
monic solution in this case is known even for quasilinear generalizations 
of fl7A0|) . see jPSOS] . 
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